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We construct minimum effective models to investigate the pairing symmetry in the newly discovered quasi- 
one-dimensional superconductor K 2 Cr 3 As 3 . We show that a minimum three-band model based on the d z 2 , d xy 
and d x 2_ y 2 orbitals of one Cr sublattice can capture the band structures near Fermi surfaces. In both weak and 
strong coupling limits, the standard random phase approximation (RPA) and mean-field solutions consistently 
yield the triplet p, -wave pairing as the leading pairing symmetry for physically realistic parameters. The triplet 
pairing is driven by the ferromagnetic fluctuations within the sublattice. The gap function of the pairing state 
possesses line gap nodes on the k z = 0 plane on the Fermi surfaces. 

PACS numbers: 74.20.Rp, 74.20.-z, 74.20.Pq 


I. INTRODUCTION: 

Searching for triplet superconductivity (SC) has been one 
of major research efforts recently partly due to its intrin¬ 
sic connection to topologically related physics and quantum 
computations. Unconventional triplet SC is normally driven 
by electron-electron interactions. Until now, the candidates 
for triplet pairing superconductors are rather limited. For 
quasi-two dimensional electron systems, Sr 2 RuC> 4 1,2 is con¬ 
sidered to be a possible good candidate. For quasi-one- 
dimensional(QlD) systems, since the electron-electron corre¬ 
lation can be further enhanced, triplet pairing has also been 
proposed in previously discovered Q1D systems, such as 
Bechgaard salts 3,4 , TUMogSeg 5 and Lio.gMogOiz 6-9 . How¬ 
ever, in these Q1D materials, the correlation may be offset 
by the fact that their electronic physics is not attributed to the 
3 d orbitals which are known to produce the strongest correla¬ 
tion effect due to small orbital size. Thus, the newly synthe¬ 
sized Q1D superconductors A 2 Cr 3 As 3 (A= K.Rb.C’s ) 10 12 be¬ 
come intriguing materials to explore possible triplet pairing 
because their electronic structures are attributed to the 3d or¬ 
bitals of chromium atoms. 

The chromium-based superconductors have attracted great 
research interests recently l0-14 . Their superconducting mech¬ 
anism shown in many experimental results 13-16 is very likely 
related to electron-electron interactions. Among the families 
of the chromium-based superconductors, the Q1D A 2 Cr 3 As 3 
(A=K,Rb,Cs) superconductors family 10-12 , which has T c up 
to 6.1 K 10 at ambient pressure, is particularly intriguing. 
The crystal of this family consists of alkali metal ions sepa¬ 
rated [(Cr 3 As 3 ) 2_ ]oo double-walled subnanotubes. The elec¬ 
tronic structures calculated from Density Functional Theory 
(DFT) suggests that the 3d orbitals of Cr make major con¬ 
tribution near Fermi surfaces(FSs) 19 - 20 . The large electronic 
specific-heat coefficient 10 , the non-fermi liquid behavior in 
the normal state 10 and the abnormality in NMR relaxation- 
rate 17 all suggest possible strong electron correlations in the 
family. The linearly temperature-dependent London pene¬ 
tration depth 18 , the lack of the Hebel-Slichter peak in NMR 
relaxation-rate upon the superconducting transition 17 , and the 


\///-dependence of C v in the mixed state 1 1 suggest unconven¬ 
tional SC with line nodes in the system. More importantly, the 
upper critical field largely exceeds the Pauli limit 10-12 , imply¬ 
ing possible triplet pairings in the family. 

The triplet pairings were also argued theoretically 20 . The 
theoretical calculations predicted that the materials are very 
close to a new in-out coplanar (IOP) magnetic state. The mag¬ 
netic state is caused by the strong antiferromagnetic couplings 
between the nearest neighbour (NN) Cr atoms and the weak 
ferromagnetic couplings between the next nearest neighbour 
(NNN) ones. Intriguingly, the magnetic fluctuations associ¬ 
ated with such a magnetic state can favor triplet pairing. A 
unit cell of A 2 Q 3 AS 3 is composed of two types of inequiva¬ 
lent Cr atoms which will be referred as Crl and Cr2 atoms 
as shown in Fig. 1, which form two sublattices, A and B, 
respectively. It is shown that Q1D electronic structures are 
dominantly contributed by only the Cr2 sublattice 20 . Within 
the Cr2 sublattice, the magnetic fluctuations are ferromag- 
netic(FM), which are caused by both the weak NNN ferro¬ 
magnetic couplings and effective ferromagnetic fluctuations 
generated by the NN antiferromagnetism. Thus the triplet 
pairing is very possible in A 2 Cr 3 As 3 . Furthermore, the ma¬ 
terials do not have space inversion symmetry but it still has 
the mirror plane symmetry with respect to the x — y plane. 
Therefore, along the Q1D chain direction, the triplet p z -wave 
pairing, which is odd with respect to the mirror plane reflec¬ 
tion, is a natural triplet candidate classified by the irreducible 
representations of the lattice symmetry. 

In this paper, we construct low energy effective models to 
study the pairing symmetry of K 2 Cr 3 As 3 . First, we construct 
the minimum tight binding model for the material. A six- 
band model is constructed with the Cr-d z 2 , d xy and d x 2_ y 2 
orbitals of both Crl and Cr2 atoms. At low energy, this 
model can be further simplified to a three-band model with 
only the (-/-orbitals of the Cr2 atoms. Starting from the three- 
band tight-binding model, we investigate the pairing symme¬ 
tries in two types of standard models, one with the standard 
electron-electron interaction terms including onsite Hubbard 
U and Hund’s coupling ./// and the other with the effective 
FM exchange couplings. We adopt combined random phase 
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approximation (RPA) and mean-field approaches for the weak 
and strong interaction limits respectively. Both approaches 
consistently yield the p z -wave triplet pairing as the leading 
pairing symmetry for physically realistic parameters satisfy¬ 
ing < r c ~ 2 (for '-jj- > r c , the RPA yields on-site inter¬ 
orbital / y 3 _ 3 T 2 y -wave pairing, driven by the Hund’s rule cou¬ 
pling), driven by intra-sublattice FM spin fluctuations. The 
gap function of this pairing possesses line gap nodes on the 
= 0 plane on the FS. This pairing is essentially Q1D and 
the pairing amplitudes on the Q1D j3 band dominates those 
on the Q1D a band and the 3D 7 - band. Finally, we shall ad¬ 
dress the experimental relevance of the p 2 -wave triplet pairing 
obtained here. 


II. MODEL AND APPROACH 


From the bandstructure of DFT calculations 19 - 20 , it’s found 
that there are five bands near the Fermi level(FL) and the orib- 
tal characters near the Fermi level are mainly of Cr- d z 2 (A\ 
symmetry) and ( d xy , d x 2_ y 2) ( E ' symmetry) components. To 
obtain an effective model, we consider a virtual atom at the 
center of each Cr triangle with three orbitals, which have the 
same symmetries as A\ and E’, as shown in Fig.l. Now 
the structure of K 2 Cr 3 As 3 is simplified as parallel placed Cr- 
chains, with each chain containing the A and B sublattices, 
which are inequivalent due to the asymmetry of K atoms and 
the absence of inversion symmetry. The point group of this 
new lattice model is the same with the original lattice and it is 
/>;/, not !){■,),, where inversion and six-fold rotational opera¬ 
tions are absent. The tight-binding Hamiltonian thus obtained 
in the momentum space can be written as, 

<B= E ^"( k ) C U-(k)Cn-(k), (1) 

k mnfiua 

where m/n = A, B labels sublattice, fjJv = 1 , • • ■ ,3 labels 
orbital and cr =f, 4 , labels spin, cj^^k) creates a spin-cr 
electron in the orbital // in the m-th sublattice with momentum 
k. The elements of the matrix /t"”' (k) and the resulting band 
structure of this model are given in Appendix.A. The band 
structure of this six-band tight-binding model well captures 
the main characters of that of DFT at low energy. 

Notice that only three of the six bands cross the Fermi sur¬ 
face, which provides possibility for us to further simplify our 
model to a three-band one. This is accomplished by hybridiz¬ 
ing each d-orbital from Crl and that from Cr2 into a lower en¬ 
ergy d-orbital, and we finally arrive at a three-orbital model. 
Neglecting spin-orbit coupling, this model is described by the 
following tight-binding Hamiltonian 



FIG. 1. (color online). The lattice model with virtual atoms. The pur¬ 
ple and green circles on each lattice site represent the virtual atoms 
of Crl and Cr2, respectively. 



FIG. 2. (color online). Band structure and Fermi surfaces, (a).The 
band structures of DFT (gray lines) and three-band tight-binding 
modelled lines). The a- and /3-bands are Q1D and the 7 -band is 
3D. (b). FSs of the tight-binding model. 


Ht B = E V(k)E(k) Ci , CT (k). (2) 


Here, the orbital index ji/v = 1, • • • ,3 represent the d z 2 for 
1, the d xy for 2, and the d x 2_ y 2 for 3, respectively. The matrix 
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TABLE I. The hopping parameters in the three-orbital tight-binding model. The onsite energies of A\ and E' orbitals are ei = 2.0939 and 
£2 = 2.0997. They are given in unit of eV. 


aB 

V 

i = z, 2 

i = z, 4 

i = z 1 6 i = ly / = 2y i = 11 y i = 22 y i = 12 y 

i = y 

i = yz 

a/3= 11 
a(3=22 
a(3=\2 

0.2841 

0.1816 

-0.0505 

-0.0469 

-0.0252 

0.0034 0.0216 -0.0135 0.0064 

0.0175 0.0095 

-0.0151 

-0.0063 



FIG. 3. (color online). Density of states and susceptibility, (a). The 
DOS (for one spin specie) of the three-band tight-binding model, 
(b). Momentum-dependence of the eigen susceptibility x (q) along 
the high-symmetry lines in the Brilloiun-Zone. The temperature is 
T = O.OOleV. 


elements h ll „{ k) are given by 

h\\ = ei + (sj 1 + 2s]^ z cos2z)(2cos2y + 4cosa;cost/) 

3 

+2^sl] 2j co S 2 jz , 
i=i 

h \2 = 2islySin2y + 2is\ y sinycosx + 2v / 3s2ySint/sina;, 
hai = h* 12 , 

hi 3 = 2slyCOs2y — 2iy/3s\ y cosysinx — 2sl y cosycosx , 

h 31 = /i* 3 , 

h 22 = e 2 + 2sll y cos2y + (sll y + 3sll y )cosxcosy 
2 

+2 S l%.i COs2 d z , 

i=1 

h -23 = V3(sjly - ^22y)sinxsin// + 2isll y sm2y, 

—4iSi2yCOSXsmy, 

^-32 = ^23> 

h 3 3 = e 2 + 2sH y cos2y + (3 sf ly + sll y )cosxcosy, 

2 

+ 2 X] s l%j cos2 j z - (3) 

2=1 

Here x = ^k x a 0 , y = \k y a a and z = \k z c 0 , where a 0 ,c 0 
are lattice constants within the xy- plane and along the z-axis 
respectively. The tight-binding parameters e/ - 2 and s“' 3 are 
shown in Table I. 

The band structure of the above introduced three-band 
tight-binding model is shown in Fig. 2(a), in comparison with 
that of DFT. The chemical potential of our tight-binding 
model is y c = 2.3105eV, leading to a band filling of 4- 


electrons per unit cell. From the comparison, it’s clear that 
the three-band tight-binding band structure well captures the 
main characters of the DFT band structure at low energy. Two 
ID bands (marked as a and /3 in Fig. 2 (a)) and one 3D band 
(marked as 7 ) are found near the FL. The two ID bands are 
degenerate along the T — A line. The FSs of the tight-binding 
model are shown in Fig. 2(b), where the FS sheets are marked 
corresponding to the three bands. The dominating orbital 
component on the a- FS is d z 2, while those on the / 3 - and 
7- FSs are (d xy ,d x 2_ y 2) (see Appendix.A). Due to the flat¬ 
ness of the two ID bands near the FS as shown in Fig. 2(b), 
large density of states (DOS) is obtained at the Fermi level, as 
shown in Fig. 3(a). 

To generate superconductivity, interactions are necessary. 
Here we focus on the possible pairing mechanism driven by 
electron-electron interactions. Since the interaction strength 
of the .‘ki-electrons of Cr might be moderate, we shall in¬ 
vestigate both the weak-coupling and strong-coupling limits 
of the electron interactions in the following, with both limits 
giving consistent results. For the weak coupling limits, we 
adopted the Hubbard-Hund’s type of on-site interactions and 
performed the standard multi-orbital RPA : 30 calculations to 
study the pairing symmetry. For the strong coupling limit, 
we adopted the super-exchange type of interactions and per¬ 
formed a mean-field analysis on the t — J type of model in¬ 
stead. 


III. THE WEAK COUPLING LIMIT 


In the weak coupling limit, the following on-site Hubbard- 
Hund’s interactions are adopted 

-^H—H U ^ ^ "F F ^ T JH ^ ^ 

i[i i,/x< v 

c !/iCT C L CT ,c w' C «!'C7- + ( c 4/4 C l/[4 C *i4 c i!'t + h.C.) . 

era' 

(4) 

Here, the [/-term, the F-term and the .///-term denote the 
intra-orbital, inter-orbital Hubbard-repulsion and the Hund’s 
rule coupling as well as the pair hopping. For simplicity, we 
have assumed that the coupling constants U, V , ./// do not de¬ 
pend on the orbital indices in our effective model. General 
symmetry argument requires U = V + 2 .///. Note that physi¬ 
cally the orbital bases d z 2,d xy and d x 2_ y 2 building our effec¬ 
tive model are not the maximumly-localized Wannier wave 
functions around a real single Cr-atom, but the mixture of the 
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three atomic wave functions around the virtual center. There¬ 
fore, the on-site interactions on these delocalized orbital bases 
are much weaker than that on a real single atom. Further¬ 
more, the interaction parameters U, V and ./// should be ex¬ 
pressed by those interaction parameters on the local Wannier 
basis, which in general can be estimated from first principle 
approaches such as cRPA, although such estimation is usually 
coarse and strongly approach-dependent. Here, we let U and 
the ratio r = Jh/U to be tuning parameters and study the 
parameter-dependence of the results. We left the determina¬ 
tion of the interaction parameters for future study. 

To study the pairing symmetry of this Hubbard- 
Hund model, we adopt the standard multi-orbital RPA 
approach 21-30 . We first define and calculate the bare suscep¬ 
tibility tensor Xi^il ’ l2 (q, t). After that, the renormalized 

spin(s) or charge(c) susceptibilities X/*^ 1 ’* 2 ( < T T ) are ob¬ 
tained in the RPA level. Then, through exchanging the spin 
or charge fluctuations, we obtain the effective pairing poten¬ 
tial V“^(k,q). Finally, solving the linearized gap equation 
for F" /3 (k. q) as an eigenvalue problem, we obtain the lead¬ 
ing pairing gap function as the eigenvector corresponding to 
the largest eigenvalue A. Generally, the RPA only works well 
for weak interactions. 

For the weak-coupling limit, we adopt multi-orbital RPA 
approach 21-30 to study the pairing symmetry of the Hubbard- 
Hund’s model of the system. The Hamiltonian of our model 
is, 


H = Htb + -Hh-h, (5) 

where the interaction term fTn-H is give by Eq.(4). 


A. Susceptibilities and magnetic phase 


Let’s define the following bare susceptibility tensor. 


xE 1 ' 2 (q.T) (k lf r)c /a (ki + q,r) 

ki,k 2 


c^(k 2 +q,0)c; 4 (k 2 ,0)) 0 , (6) 

with Zj/j—i ... 4 ) denoting orbital indices. The explicit formal¬ 
ism of is given by 

xJS’* 2 (q- iw n) = ir 5Z O k )C*( k )€( k + < l) 


k,o;,/3 


^(k + q) 


n F(e k +q ) - M £ k) 


,(7) 


lUJ n 


- e k+q 


where a//3 = 1,..., 3 are band indices, eg and (k) are the 
a —th eigenvalue and eigenvector of the h{ k) matrix respec¬ 
tively and riF is the Fermi-Dirac distribution function. 

To detect the pattern of the dominating spin correlation 
in the system, we define the static susceptibility matrix 
Xq {l,ni) = Xm]m(q,iu} n = 0). The largest eigenvalue of this 
Hermitian matrix for each momentum q is defined as x(q), 
which represents the eigen susceptibility in the strongest chan¬ 


TABLE II. The ten possible pairing symmetries for K 2 CX 3 AS 3 , 
among which five are spin-singlet while the left are spin-triplet. 


singlet 

triplet 

s 

Pz 

(d x 2 — y2 ) d X y^) ( d x 

2 y 2 7 d* X y ) * Pz 

( Px,Py) -Pz 

(Px,Py) 

fx 3 —3xy 2 ‘ Pz 

fx 3 —3xy 2 

fy 3 —3x 2 y ‘ Pz 

fy 3 —3x 2 y 


nel. The distribution of x(q) in the Brillouin-Zone reflects 
the pattern of spin correlations in the system. Fig.3(b) shows 
the momentum-dependence of x(q) along the high-symmetry 
lines, where the peak at the F-point reveals strong FM spin 
correlations in the model. Note that in the building of our 
model, only one sublattice is taken from real material, which 
thus means that the spin correlation pattern found here is intra¬ 
sublattice FM. 

When interactions turn on, the renormalized spin(s) or 
charge(c) susceptibilities are obtained in the RPA level as 


X (s ’ c) (q, iv) = [/ =F X (0) (q, iv) U^] X (0) (q, iv) (8) 


Here x^’^ (q, iv n ), (q, w») and U( s,c ' ) are operated as 
9x9 matrices (the upper or lower two indices are viewed 
as one number) and the matrix U^ s,c ' ) are given in the 
Appendix.B Clearly, the repulsive Hubbard-interactions en- 
hance(suppress) spin(charge) susceptibility. Note that the 
RPA approach only works for weak U ^ when all the eigen¬ 
values of the denominator matrix /—X® (q, iv) U ^ are pos¬ 
itive. At the critical interaction strength, the lowest eigenvalue 
of that matrix touches zero, which leads to divergence of the 
spin susceptibility and invalidates the RPA, suggesting the for¬ 
mation of magnetic order. The critical strength U c as function 
of Jh/U is given in the phase-diagram in Fig. 6 as the phase 
boundary separating the superconducting phase and the mag¬ 
netic phase. For U > U c , intra-sublattice FM phase emerges 
since the bare susceptibility shown in Fig.3(b) peaks at the F- 
point (the renormalized susceptibility also peaks there), con¬ 
sistent with the “IOP” phase obtained in Ref.20. 


B. Pairing symmetry study 

When the interaction strength U < U c , there can be 
short-ranged spin or charge fluctuations in the system. Let’s 
consider a Cooper pair with momentum/orbital (k 't, —k's), 
which could be scattered to (kp, —kg) via exchanging spin 
or charge fluctuations. This second-order perturbation pro¬ 
cess will contribute an effective interaction vertex F 2 ) 2 (k. k'). 
Considering only intra-band pairings, we project the effective 
interaction F^(k,k') onto the LSs and obtain an effective 
pairing potential H“^(k, k'), where a/f3 = 1, • • • , 3 are band 
indices. The explicit formula of T^k, k') and l /0:/:<1 (k. k') 
are given in the Appendix.B. 
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Finally, one should solve the following linearized gap equa¬ 
tion to find the pairing symmetry. 


1 

( 2^)3 



V a/ 3 (k, k') 

< 4 ( k ') 


A^(k') 


AA a (k). 


(9) 

Here, the integration and summation are along various FS 
patches. The Vp( k') is Fermi velocity and k'| represents the 
component within patch /3. This equation can be looked upon 
as an eigenvalue problem, where the eigenvector A a (k) rep¬ 
resents the relative gap function near T c , and the eigenvalue 
A determines T c via T c =cut off energy e -1 / A . The leading 
pairing symmetry is determined by the largest eigenvalue A. 
Generally, the RPA study on pairing symmetry only works 
well for U « U c , which will be our focus in the following. 


The eigenvector(s) A,, (k) for each eigenvalue A introduced 
above form an irreducible representation of the D 3 ^ point 
group of K 2 CT 3 AS 3 . Altogether, there are 10 different pair¬ 
ing symmetries, including 5 singlet pairings and 5 triplet 
ones, which are listed in Table.II. Here, the (d x 2_ y 2,d xy ), 
( d x 2_ y 2,d xy ) -p z , ( p xi p y ) and (p x ,p y ) ■ p z symmetries each 
form a 2D representation of the point group with degenerate 
pairing eigenvalues. Generally, below T c , the two degener¬ 
ate eigenvectors of these representations would be mixed as 
d x 2_ y 2±id xy ,{d x 2_ y 2±id xy )-p z ,p x ±ipy and (p x ±ip y )-p z . 
There are two different /-wave symmetries, i.e. f x 3_ 3xy 2 and 
fy 3 - 3 x 2 y While the gap functions of both /— waves change 
signs with every 60° degree rotation, their nodal lines are dif¬ 
ferent by 30° rotation. Due to the absence of spin orbital cou¬ 
pling (SOC) in our model, each of the five triplet pairings pos¬ 
sesses three degenerate components, i.e. ft> (tl + lt)> 
which we shall not distinguish in the following unless oth¬ 
erwise stated. Note that the concrete formalism of the gap 
function of each symmetry is not given a prior to the RPA cal¬ 
culation. Instead, they should be solved from the gap equation 
(9). 


Our RPA calculations were performed on 18 x 18 x 200 
and 30 x 30 x 300 lattices with periodic boundary condi¬ 
tion. The consistence between the two results indicates our 
RPA result converges in the thermodynamic limit. The 17- 
dependence of the largest eigenvalue A for the five stronger 
pairing symmetries is shown in Fig. 4 for (a). Jh = 0.217 
and (b). Jh = 0.417. Clearly, all these eigenvalues increase 
promptly with the enhancement of U and would diverge for 
U —> U c . For Jh = 0.217 shown in (a), the p z -wave is the 
leading pairing symmetry and dominates other symmetries; 
for Jh = 0.417 shown in (b), the /j / 3 _ 3 x 2 y -wave becomes the 
leading one, with the (p x ,p y )- and p z - waves to be close can¬ 
didates. Fig. 5(a) shows the ^-dependence of these eigen¬ 
values for fixed U = O.leV. Obviously, the spin-triplet p z , 
f y 3 - 3 X 2 y and ( p x ,Py ) symmetries dominate other symmetries 
in the whole range of r = While for r < r c (~ |) the 
p z -wave symmetry dominates other symmetries, for r > r c 
the / y 3_ 3 a ,2 y -wave symmetry becomes the leading symmetry 
with the ( PxiPy ) and p 2 -wave states to be close candidates. 




U(eV) ' ' U(eV) 

FIG. 4. (color online). The 17-dependence of the largest eigen values 
A for different pairing symmetries for (sl).Jh = 0.21/ and (b). Jh = 
0.417. Only the five stronger symmetries are shown. 



FIG. 5. (color online). The -dependence of the largest eigen val¬ 
ues A for the five stronger pairing symmetries for fixed U = O.leV 
for (a). RPA and (b). mean-field results. 


C. Triplet pairings 


It’s interesting to find from the above shown RPA re¬ 
sults that all the dominating pairing symmetries in differ¬ 
ent regimes, i.e. the p z , f y 3 _ 3 x 2 y and ( p x ,Py ) symmetries 
are triplet pairings. The ground state phase-diagram in the 
U — j- plane is shown in Fig. 6 , where three possible phases 
are present. For U > U c (which is around 0.16-0.24 eV 
and is J n /(/-dependent), the intra-sublattice FM SDW order 
emerges. For U < U c , two triplet pairings, i.e. the p z and 
f y 3 _ 3 x 2 y -waves emerge, which are separated by the critical 
value r c (which is 17-dependent) for r = jj-. For r below 
or above r c , the p z - or f y 3 -wave pairing is the leading 
pairing symmetry respectively. In the limit of U —» 0, we have 
r c —> |; with the enhancement of U, r c first decreases slightly 
and then drops near the SDW critical point. In our focused 
regime of 17 << 17 c in RPA, we have r c » While both 
gap functions have node lines, they are along k y = ±\/3 k x , 0 
directions for the f y 3 _ 3 x 2 y - wave pairing and in the k z — 0 
plane for the p z -wave pairing. 

Note that the dominating part of the leading f y 3_ 3x 2 y - and 
sub-leading (p x ,p y )- wave pairings in the regime r > r c « | 
are on-site inter-orbital pairings driven by the Hund’s rule cou¬ 
pling. To clarify this point, let’s turn off the second-order per¬ 
turbation process of RPA, and study the problem in the mean- 
field level. The Hubbard-Hund’s interaction (4) can be decou- 
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FIG. 6 . (color online). The ground state phase-diagram in the U— 7/ 
plane. 


pled in the real space as the sum of an interorbital pairing part 

-Winter = (V + Jh) ^ W EE W 

i,fi< v 

+ (V — ./#) ^ ' ^lii/crcT 00 ‘E^o- (*) 

i,fi<.u,a 

+ E EE«EE«] no) 

and an intraorbital pairing part 31 

#intra = £7 E/ E/d 7 W EmIT (*) 
i[i 

+ jh y At (o a,*, w, (ii) 


whete ^uvaa' (^) — CificrCiva' and — 

^ ^A^-j-p ± A^^. Clearly, for repulsive interactions, 
energetically favored mean-field pairings are only possible 
in the regime Jh > V (or equivalently r = J T j > r c = '!, 
upon U = V + 2 Jh)- The possible pairings channels include 
two on-site inter-orbital triplet pairings, which are the leading 
/y3-3x^-wave: A 23)((7(7)T4-+4 , t ) and the subleading (p x ,p y )- 

doublets: (A 12)(<7<Tin+4 . t) , A 13i(crcrjn+4 . t) ), whose pairing 

eigenvalues as functions of are shown in Fig.5(b). The 
real-space pairing pattern of the mean-field f y 3_ 3x 2 y -wave 
pairing is shown in Fig.7(a) (only the // pairing component 
is shown), where the pairing is between the d xy and d x 2_ y 2 
orbitals. Comparing Fig.5(a) and (b), the FM spin fluctuations 
considered in RPA further enhance the two triplet pairings 
and slightly reduce r c , besides extending the pairing coherent 
length. 

The distribution of the relative gap function of the leading 
/-wave pairing obtained by RPA for U = O.leV, Jh = 0.4(7 
is shown on the FSs in Fig.8(a),(b) and (c). This orbital-singlet 
spin-triplet pairing state is mirror-reflection even about the 
k z = 0 plane, and its gap function changes sign over every 


60° degree rotation about the 2 -axis, causing line gap nodes 
in the direction of k y = 0, ±\/3 k x . While the pairing ampli¬ 
tudes on the /3- and 7 -FSs are comparable, that on the cv-FS 
is much weaker as the dominating orbital component on that 
FS is d z 2 , which nearly doesn’t take part in the pairing, as 
suggested by the above mean-field analysis. 

In real material, the situation ^ > r c ~ | can hardly 
be realized because the orbitals in our effective models are 
the molecule orbitals that are equally distributed on the three 
Cr2 atoms and thus the effective value of Jh should be small. 
Therefore, the realized pairing symmetry in fCCr-iAsa might 
more probably be the p z -wave pairing displayed in the regime 
r < r c in the phase-diagram Fig. 6 , which is consistent with 
the mean-field study below on the t—J model of the system. 
Since the Hund’s rule coupling driven pairings don’t include 
the p z -wave, as suggested by Fig.5(b), this pairing is purely 
driven by the intra-sublattice FM spin fluctuations reflected in 
Fig.3(b). 

The distribution of the relative gap function of the leading 
p z -wave pairing obtained by RPA for U = O.leV, Jh = 0.1(7 
is shown on the FSs in Fig. 8 (d), (e) and (f). This triplet pair¬ 
ing is mirror-reflection odd about the k z = 0 plane, caus¬ 
ing line gap nodes on that plane on the 3D 7 -FS. The gap 
function of this pairing is isotropic on the FS sheets with a 
fixed k z and its k x - and k v -dependence is weak; particularly 
it doesn’t change sign all over the regime k z > 0 (or k z < 0 ) 
on the FSs, causing no extra gap nodes. Obviously, the gap 
amplitude on the Q1D ,3-FS dominates those of a- and 7 - 
FSs, suggesting the Q1D nature of the obtained p z -wave pair¬ 
ing. Considering the orbital components of /3-FS, this pairing 
takes place mainly within the d xy and d x 2_ y 2 orbitals. To 
exhibit the fc- -dcpcndcnce of the obtained A a (k), we plot in 
Fig.7(b) (A a (k)) FS .i.e. the averaged gap function on the FS 
as function of the fixed k z , in comparison with that obtained 
in the following strong-coupling limit in Fig.7(c). The consis¬ 
tency between the results obtained from the two approaches 
are quite good. 


IV. THE STRONG COUPLING LIMIT 


In the strong coupling limit, we consider super-exchange 
type of interactions. From DFT calculations on K 2 Cr 3 As 3 20 , 
the NN and NNN exchange couplings along the Cr-chain are 
antiferromagnetic (AFM) and FM respectively. As only the 
Cr2 atoms are included in our three-band model, the super¬ 
exchange interactions between adjacent Cr2 along the chain 
direction (i.e. the 2 -direction) are actually FM, which pro¬ 
vides us the following exchange interaction term 

Hj = -J 2 y S ai -S aj , (12) 

with J 2 > 0. Here a £ {d z 2 ,d xy ,d x 2 _ y 2 j, ( ij) denotes 
NN bond along the 2 -direction in our model but NNN bond 
in real materials and S Qi = \ c aia tT ™’ c cdcr' is the lo¬ 
cal spin operator. The inter-chain exchange coupling has been 
neglected here due to the Q1D character of the material. Note 
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FIG. 7. (color online), (a). Dominating pairing pattern in the real space for the / 3_ 3a .2 y -wave pairing shown in Fig.(6). Note that only the 
pairing component is shown. The k~ -dependence of the relative gap function of p z -wave pairing averaged on the FSs for (b). RPA result for 
the Hubbard-Hund model with Jh/U = 0.1 and (c) Mean-field result for the t—J model. 


that only intra-orbital FM exchange coupling is considered be¬ 
cause of the “IOP” spin configuration introduced in Ref.20. 
Adding this interaction term into the tight-binding model, we 
obtain an effective t—J model. Since the band filling here is 
far from making the system a Mott-insulator even in the large 
U limit, we shall omit the no-double-occupance constraint and 
perform a mean-field analysis on the model. 

We consider the super-exchange interaction Eq.(12). The 
Hamiltonian of our model is, 


H = H tb + Hj. (13) 


Since this material is far from Mott-insulator, we neglect the 
no-double-occupance constraint on this t — J model and per¬ 
form a mean-field study on it. 

In terms of real-space pairing operators, this interaction 
Hamiltonian can be written as. 


— 4^2 EE 4 ^2 

{ij) 


T 1 T 


EE 1 ijn,'Y- L n» m 
(D').m t 


pi =-U 

ij ij. 


V2 K 


,t t 
"*/4 j/4 


J J 
ipt 1/4 


J J ) 

Vl 1/4'’ 


T' = ! J-(P P +r f P ) 

iji^n \ I */4 1/4 ~ Vl i/*t' 

P J 


C */4 C l/4 


7 = 1 . 

7 = 0 , 

7 = - 1 . 


(14) 


Here P, :ill is the spin-singlet pairing operator and is the 

spin-triplet one with S z = 7 . From Eq. (14), it’s clear that for 
J 2 > 0 (J 2 < 0), spin-triplet (singlet) pairing will be favor¬ 
able. In K 2 Cr 3 As 3 , FM intra-sublattice spin correlation leads 
to J 2 > 0, thus we only consider the triplet pairing terms. 
Due to the absence of SOC in our model, the pairing states 
in S z = —1,0,1 channels are degenerate, which allows us to 
only consider the S z = 0 channel in the following. 

Defining A M = -f J 2 (T i>i+5iMi o)/2\/2 as the (real) pairing 
order parameters, the above Hj term can be decoupled as 


H j - EE {~ iA A,o + h - c ) + EE |a m | 2 - 

(»i>.M 2 v 


(15) 

In combination with the tight-binding part, the total mean- 


field Hamiltonian in the momentum space is 

p j\r 

H mi = £ T-t(k)7T(k)vk(k) + — EE | A 


2 

Ml ’ 


H(k) - ( ¥ k) ^ (k) ^ 

( } VA|,(k) 


Atl(k) = 


Ai(k) 


A 2 (k) 


As(k) 


(16) 


where (k) = [c{ f (k), 4 t (k), 4 f (k), c u (-k), c 2i (-k), 
c 3 f( — k)] and A^(k) = A^sin k z . Note that a constant 
Sk//[^W(k) — Pc], which has no contribution to the dynam¬ 
ics, has been neglected on the above. 

The above mean-field Hamiltonian can be solved 
via diagonalizing // (k) by a unitary transformation: 
£/l(k)fT(k)(/(k) = Diag(i?i), with the eigenvalues 

Ei = —Ei + 3 > 0 (i = 1,2,3). The self-consistent gap 
equations are. 


a . = ^EE ^+3, a (k)^,a(k)/[^(k)]sinfe, 

k a—1 

(17) 

where f{E) is the Fermi-Dirac distribution function and we 
only study the case of zero temperature. Solving these gap 
equations, we plot the ground state energy Eg and the three 
gap amplitudes A (J (// = l,--- ,3) as function of J 2 in 
Fig. 9(a) and (b) respectively. Clearly, the gap amplitudes of 
d xy and d x 2_ y 2 orbitals, belonging to the E' irreducible rep¬ 
resentation, are equal and dominate that of d z 2, which is con¬ 
sistent with the RPA result. Furthermore, there is no phase 
difference between the SC order parameters of the three or¬ 
bitals. 

Fig. 10 shows the distribution of the pairing gap function 
•^2 (Cc4 (—k) c at (k)) = J 2 X]^(c^(-k)c pt (k))|C At a(k )| 2 = 
4^ A A1 |^ AtQ (k)| 2 sin(/c 2 ) overthe FSs in the Brillioun-Zone. 
Here a is the band index and only intra-band pairing is shown. 
Clearly, due to the same sign among different A ;l , there is no 
sign change of the gap function all over the regime k z > 0 
and gap nodes only exist in the k z = 0 plane in the 3D 7 - 
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FIG. 8. (color online). The relative gap function distribution over the three FSs for the f y 3_ 3x 2 y wave state(top panel) and the p z wave 
state(bottom panel) in RPA calculations. Fig. (a) (d), (b) (e) and (c) (f) are for a, ft and 7 FSs. With U = 0.1 eV, the adopted Hund couplings 
are Jh = 0.417 and Jh = 0.117 for f y 3_ 3x 2 y and p, wave states in the calculations. The gap functions are plotted on the 30x30x300 
k-mesh. Note that the strong local oscillations of the gap are artifact of plotting due to finite lattice. 



that of E'{d xy , d^.^-orbitals. Due to the phase difference 
and orbital characters, the gaps on the a- and 7- FSs are 
strongly anisotropic and line nodes appear on the a- FS; but 
the gap on the /3- FS is quite uniform. Compared with the for¬ 
mer case, there is an additional coherence peak in DOS. This 
state, however, is energetically unfavorable in our results, due 
to the extra gap nodes on the FS. 32 . 


V. DISCUSSION AND CONCLUSION 


FIG. 9. (color online). The ground state energy Eg (a) and gap 
amplitudes A M (b) for triplet pairing in S z = 0 channel as a function 
of,/ 2 . 


band. Obviously, the pairing symmetry of this superconduct¬ 
ing state is p z since its gap function is reflection odd respect 
to the mirror plane k z = 0 and is rotation symmetric about 
the 2 -axis. What’s more, due to the dominance of A 2i 3 over 
Ai, the pairing amplitude on the Q1D /3-band is obviously 
stronger than that on the o-hand. Fig. 11(a) shows the gap 
distribution over the FSs in the k x — k z plane with k y = 0 
and Fig.7(c) shows the averaged gap function (A a (k ))pg on 
the FSs with fixed k z in comparison with the same quantity 
obtained from RPA in Fig.7(b). The obtained gap structure 
here is well consistent with the results of RPA introduced 
above. In Fig. 11(b), we plot the DOS in the superconduct¬ 
ing state with Ai = 2.5meV, A 2 = A3 = lOmeV obtained 
for J 2 = 0.4eV, scaled by a factor of 2 for clarity. The DOS 
displays a V-shaped structure, which can be detected by tun¬ 
neling experiments. 

Note that there can be another p 2 -wave solution of the gap 
equation with the gap sign of ,4' — d z 2 -orbital different from 


In conclusion, we have performed combined weak- 
coupling RPA and strong coupling mean-field study on the 
pairing symmetry of K 2 Cr 3 As 3 , starting from the three-band 
tight-binding model. The weak-coupling RPA study on the 
Hubbard-Hund model reveals two possible triplet pairings in 
different parameter regimes. For r = 7^ > r c « on¬ 
site inter-orbital pairing (between d xy and d x 2_ y 2) triplet SC 
would be driven by the Hund’s mle coupling. This pairing 
has nodal lines in the direction of k y = 0, ±\/3 k x . However, 
as the orbitals in our effective models are the molecule or¬ 
bitals that are equally distributed on the three Cr2 atoms, the 
effective value of Jh should be small. Therefore, r < r c 
is a more realistical parameter. For r < r c , the Q1D p z - 
wave pairing would be driven by the strong intra-sublattice 
FM spin fluctuations. The p z - wave pairing state obtained by 
RPA for the Hubbard-Hund model is in well consistency with 
that obtained by mean-field study on the t — J model in the 
strong coupling limit. The gap function of this pairing state 
doesn’t change sign in the k z > 0 (or k z < 0) regime on the 
FSs and the only line gap nodes lie on the k z = 0 plane on the 
3D 7 -FS. The pairing amplitude on the Q1D /?- FS dominates 
those on the Q1D a- FS and the 3D 7 - FS, revealing the dom¬ 
inating part of this pairing to be the Q1D pairing within the 
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FIG. 10. (color online). The gap function distribution over the three FSs with A^| = 0.0025eV, A^_ a = A 44 = O.OleV, which are twice 
the gap values obtained for J 2 = 0.4eV. Fig. (a), (b) and (c) are for a, /? and 7 FSs. 



FIG. 11. (color online). The SC gap function distribution over the 
FSs in the k y = 0 plane (a) and density of states (for one spin specie) 
(b) with A^i = 0.0025eV, AT£ 2 = A 44 = O.OleV, which are 
twice the gap values obtained for J 2 = 0.4eV. The red and blue lines 
represent DOS for normal and superconducting states, respectively. 



J„/U 

FIG. 12. (color online). The ^--dependence of the largest eigen 
values A for the p z - and / 3_ 3x 2 y - wave pairings for fixed U = 0.1 
eV and X ao = 10 meV for the S z = 0 (T|) component and ,31 = ±1 
one. 


E'{d X y, d x 2 _ y 2 ) orbitals along the z-axis. 

In K 2 Cr 3 As 3 , the weak SOC will lift up the degeneracy 
among the three components of the triplet pairings. The only 
symmetry allowed on-site SOC term added to our three-band 
tight-binding model takes the following formalism. 


Hsoc = iKo ^ ° C L( k ) c 3<r( k ) - 4<r( k ) c 2 CT ( k ) 


, (18) 


kcr 


where a coupling constant A = lOmeV is adopted in our cal¬ 
culations. Classified according to the S z quantum number 
of Cooper pairing, possible pairing components include the 
S z = ±1 (i.e. j"MI) components and the S z = 0 (i.e. fl) 
component (note that due to the absence of inversion symme¬ 
try in K 2 Cr 3 As 3 , the triplet T| + IT component in principle 
is allowed to be mixed with the singlet f| — || component 33 , 
although such mixing is weak due to the weak SOC in 3d elec¬ 
tronic materials). The iff- dependence of the largest pairing 
eigenvalues A of the S z = ±1 and the S z = 0 components of 
the p z - and the f y 3_ 3x 2 y - wave pairings obtained from RPA 
(details of the formula of RPA with SOC can be found in Ap- 
pendix.B.) is shown in Fig. 12, with fixed U = O.leV. Clearly, 
for both symmetries, the S z = 0 component slightly wins over 
the S z = ±1 ones. Our mean-field calculations for the t — J 


model also energetically prefer the S z — 0 component to be 
the true ground state with a small energy gain. 

Note that the gap amplitude on the Q1D /3-band dominates 
over that on the other two bands in the p=wave pairing state 
obtained here, which suggests that Cr2 is more important than 
Crl in determining the pairing properties of the system. This 
self-consistently justifies our three-band model as a good start 
point for the study of pairing symmetry. Nevertheless, to ob¬ 
tain more detailed physical properties of the system, particu¬ 
larly the magnetic property, the six-band model with both Crl 
and Cr2 involved is necessary, as both of them participate in 
the magnetization and the relative orientation between their 
magnetic moments are opposite 20 . We leave this problem for 
future study. 

On the experimental aspect, the line gap nodes of this p z - 
wave pairing manifest themselves in the specific heat 1011 , the 
NMR relaxation-rate 17 , and the penetration-depth 13 . Besides, 
a direct ARPES observation will reveal the existence of line 
gap nodes on the k z = 0 plane on the FS. The feature of triplet 
pairing in this material is consistent with the strong upper crit¬ 
ical field largely beyond the Pauli limit 10,1 1 . In addition, it will 
be reflected in the lack of obvious change in the NMR knight- 
shift upon the superconducting transition 34 . Furthermore, the 
spin triplet pairing state here should be very sensitive to non- 
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magnetic impurities. The Q1D pairing nature of this p z -wave 
SC can be detected by anisotropic coherent length. In addi¬ 
tion, a direct ARPES observation can reveal the dominance 
of the gap amplitudes on the Q1D /3-band over the other two 
bands. 
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Note added'. After the completion of the present work, we 


noticed the work by Yi Zhou, et al 35 , studying topics that 
partly overlap with the present work, but their conclusions are 
different from ours. After comparing their work with ours, 
we notice the following difference. Firstly, the band structure 
adopted in our study is different from theirs, particularly in 
the aspect of orbital components on each FS. Secondly, in the 
RPA study, we obtained the relative gap function for each pair¬ 
ing symmetry via solving the linearized gap equation, which 
guarantees that the leading gap form-factor is not missed. In 
addition, the issue of critical interaction strength in RPA has 
been thoroughly investigated in our study, which guarantees 
that the interaction parameters adopted in our calculation is 
below the critical one and the RPA is valid. 


* yangfan_blg@bit.edu.cn 
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Appendix A: The tight-binding model 


From the bandstructure of DFT, we find that there are five bands near the Fermi level and the oribtal characters are mainly d z 2, 
d xy and d x 2 _ y 2 . To obtain an effective model, we consider a virtual atom at the center of each Cr triangle with three orbitals, 
which have the same symmetry as A[ (d z 2 ), and E' ( d xy , d x 2 _ y 2 ), as shown in Fig.l in the main text. Now the structure of 
K 2 Cr 3 As 3 is simplified and every chain has two sublattices. Due to the asymmetry of K atoms and the absence of inversion 
symmetry, there are two inequivalent sublattices which are labeled as A and B respectively. The tight binding Hamiltonian is 


K k) 


/ h AA h AA h AA h AB h AB h AB \ 


°12 ,l 13 l0 U '*12 

uAA uAA uAB uAB 
/1199 f c<23 ' f'D 1 it" 


13 

h AB 

,l 23 


21 '"22 
uAA uAB uAB uAB 
11 33 ,l 31 ,l 32 ,l 33 

h BB h BB h BB 

^22 B h 23 


h 33 J 


(Al) 


where the unshown matrix elements can be obtained from the shown ones by Hermicity of the /i(k) matrix. On the above, the 
orbital indices 1, 2 and 3 denote d, 2 , d xy and d x 2_ y 2 respectively. 


Let x = ^k x ao, y = \k y ao and z = \k z Co, the Hamiltonian matrix elements are as following. 


(1). The hoppings between A\ (d z 2 ) orbital and A[(d z 2 ) orbital lead to the following matrix elements of h{ k) 


h AA/BB = ei/3 + 2^2 sl] 2i cos2iz + (s 


2s y 1 z cos2z)(2cos2y + 4cos.tcos y) 


h AB = 2 ^2 s ^2i-i cos [(2* — 1 ) 2 ] + 2s y 1 z cosz(2cos2y + 4cosa;cos y) 


i =1 

uBA _ t AB,* _ uAB 

'*11 ~ '*11 ~ '*11 • 


(2). The hoppings between A! 1 [d z 2 ) orbital and E'(d xy , d x 2 _ y 2 ) orbitals lead to 

hf^ BB = 2is{ y sm2y + 2is\ y smycosx + 2v / 3s2ySinysinx 


h. 


AA/BB _ yAA/BB,* 


hf^ BB = 2s2yCOs2y — 2iv / 3s^cost/sin.T — 2sl y cosycosx 

, AA/BB _ , AA/BB,* 

"31 ~ "'13 

h AB = 2cosz(2is\ yz sm2y + 2is\^ z sinycosx + 2v / 3s2y Z sinysinx) 

uAB _ > AB,* ,BA __ h AB j BA _ ,AB,* 

"21 — "l2 7 "12 — "12 7 '*21 — "l 


l 12 


l l2 


h AB = 2cosz(2sl yz cos2y — 2v / 3is^ 2 cosysin:r — 2s^ z cosycosx) 


’2 ~~ ""■’l yz 

uAB _ 1 AB,* , BA _ 7 AB 7 BA _ 7 AB,* 
"31 '"13 7 "13 — "13 7 "31 — " 1 


l 13 


‘13 


(A2) 


(A3) 
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TABLE III. The hopping parameters (in unit of eV) along c axis to fit the DFT results in the six-orbital model. The onsite energies are 


ei = 1.9080eV, 

e 2 = 2.0407eV, £ 3 = 1.9722eV and £ 4 

= 1.9412eV. 





a/3 

V 

i = z,l i = z, 2 

i = z, 3 

i = z, 4 

i = z, 5 

i = z, 6 

i = z, 7 

a/3=U 

a(3=22 

0.1749 0.1568 

-0.0041 0.1734 

0.0301 

0.0019 

-0.0532 

-0.0452 

-0.0097 

0 

-0.0189 

0 

0 

-0.004 


TABLE IV. The inplane hopping parameters (in unit of eV) to fit the DFT results in the six-orbital model. 

olB 

V 

i = 1 y 

i = 2 y 

i = 11 y 

i = 22 y 

i = 12 y i = y 

i = yz 

af3=\2,s\ 2 
a/3=12,sj z 
a/3=22,sf 2 
a/3=22,s 22 
a/3=ll,s 2-1 

0.0206 

0.0066 

0 

0.0133 

0.0460 

0.0061 

-0.0191 

-0.0109 

0.0062 

0.0159 

-0.0218 

-0.0031 


(3). The hoppings between E'{d xy , d x 2 _ y 2 ) orbitals and E\d xy , d x 2 _ y 2 ) orbitals lead to 

2 

h^2^ BB = £2/4 + 2 siij,cos 2 y + (s 2 f y + 3sl 2 y )cosxcosy + 2 ^ s1 2 2i cos2iz 


i=1 


h 23 BB = V3(sll y — S 22 y) s i n ‘ TS i n 2 / + %is^yS\a2y — 4is 2 2 y cosxsiny 


h 


AA/BB 

32 


= h 


AA/BB,* 

23 


h 3 ^ /BB = e 2 /4 + 2sll y cos2y + (3 sfl y + sll y )cosxcosy + 2 ^ s 22 2i cos2iz 

2 = 1 

3 

h- 2 B = 2cosz [2sfl yz cos2y + (sfj yz + 3sll yz )cosxcosy\ + 2 ^ s 22 2l _ 1 cos(2i - 1 )z 


2=1 


lBA 

l 22 

lAB 

l 23 

lAB 

l 32 

AB 

33 

BA 

33 


_ i AB,* _ i 
— ri 22 — ,l 22 


AB 


= 2cos z V3(sii yz — 522 yz )sinxsiny + 2is^ yz svn-2y — 4isi 2 yz cosxsiny 

_ 7 AB,* 7 BA _ lAB h BA _ h AB,* 

~ n 23 J rb 23 — ,l 23 i ,l 32 ~ ll 23 

3 

= 2cosz [2sf 2yz cos2y + (3 sf lyz + s 2 2 l yz )cosxcosy] + 2 ^ sl 2 2i _ 1 cos{2i - 1 )z 


i=1 


_ h AB,* _ ,AB 
— '‘33 — “33 


(A4) 


The hopping parameters are obtained by least-square-root fitting of the above tight binding model to the DFT band structure. 
These parameters are given in Table III and IV and the band structures of tight binding model and DFT are shown in Fig. 13. 
Near the Fermi level, the tight binding band fits well with that of DFT. 

As the states near the Fermi level are dominantly contributed by Cr2, we can drop Crl and further simplify our model to a 
three-band one. The Hamiltonian of the resulting tight-binding model is , 


H k) = 


/in h\2 hi3 
/*22 /*23 
/‘33 


(A5) 


Again, the unshown matrix elements are obtained by Hermicity of the h(k) matrix. The formula of the matrix elements of h{ k) 
are given in Eq.(3), with the hopping parameters in this model given in Table 1 in the main text. The band structure of this model 
is shown in Fig. 2 in the main text, in comparison with that of DFT. This tight binding band also fits well with that of DFT near 
the Fermi level. Furthermore, the orbital characters, shown in Fig. 14(a), are consistent with those of DFT 20 . Fig. 14 shows the 
weight distribution of E' orbitals on Fermi surfaces in k y = 0 plane. 
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FIG. 13. (color online). The band structures of DFT (gray lines) and tight binding(red lines). 




1 


0.8 

0.6 

0.4 



FIG. 14. (color online). Orbital characters, (a) The orbital characters of the band structure of three-band tight binding model, (b) The weight 
distribution of E'(d xy , d x 2 _ y 2 ) orbitals on Fermi surfaces in k y = 0 plane. 


Appendix B: The multi-orbital RPA approach 


The Hamitonian adopted in our calculations is 


H = H t - b + H, 


H int = U 


ifi 


int 


+ y ni^riiv + Jh c tfj.a c ti/a ,c, ‘ H<T ,c iv<7 + (c^- ^ c tn\. c ivi c iv\ + h.C.) 






(Bl) 
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FIG. 15. (color online). Feyman’s diagram for the renormalized susceptibilities in the RPA level. 


Let’s first study the case without SOC, then consider the inclusion of on-site SOC for K 2 Q 3 AS 3 . 


1. Without SOC 


Let’s define the following bare susceptibility for the non-interacting case (U = V = Jh = 0), 


jo )h,l 2 
x h,U 


(q, r ) 


^ X ( T -r c L(ki>r)ci 2 (ki+q,r)c^(k 2 +q,0)cj 4 (k 2 ,0)^ , 

V ki,k 2 


where Z, (i = 1, • • • , 4) denote orbital indices. The explicit formalism of y ,(l : ' in the momentum-frequency space is. 


(0)h,l 2 


(q, W n ) 


b S « 


(k)C 


*(k)£flk + q)£f’*(k + q) 


n F (e? 


k+q 


) - M £ k) 


k,a,/3 


I 

+ £ k+q 


(B2) 


(B3) 


where a/(3 = 1,..., 3 are band indices, eg and (k) are the a —th eigenvalue and eigenvector of the h( k) matrix respectively 
and n F is the Fermi-Dirac distribution function. 
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FIG. 16. The three processes which contribute the renormalized effective vertex considered in the RPA, with (a) the bare interation vertex and 
(b) (c) the two second order perturbation processes during which spin or charge fluctuations are exchanged between a cooper pair. 


When interactions turn on, we define the spin (x^) and charge (x <ri ) susceptibility as follow, 

(q,r) = E Y, (rr^ 1 , ffl (k 1 ,r)a 2 , ffl (k 1 +q ) r)C^ iff2 (k 2 + q ) 0)<7 j4 , O5l (k 2 ,0)), 

ki,ka,<Ti,o-2 


Xj 3 ,l4 


x[ S J: M (q,r) = E aia2 ( T ^ C '/i,<T 1 (ki ! T)C'; 2i(Tl (k 1 +q,r)C'+ CT2 (k 2 +q,0)C' ;4i<72 (k 2 ,0)J), 


kl,k2,<Tl,CT2 


*h,h ^ h ’ h (^ r ) = i E (rrC'/ it (k 1 ,r)C' / 2 i (k 1 +q,r)C+ 4 (k 2 +q, 0 )C'i 4 t (k 2 , 0 )\ , 


ki,k 2 


x\:2 )hM (q. T ) = ^ E { T r C ti^r)C h t(k x + q, r)C+ t (k 2 + q, 0 )C) a (k 2 , 0 )) . (B4) 

ki.k 2 

Note that in non-magnetic state we have x^^ = X^ s+ ^ = X^ s — X^> and when U = V = Jh = 0 we have x^ = 
X <,s> = x <{>> ■ In the RPA level, the renormalized spin/charge susceptibilities for the system are, 

-1 


X (s) (q, iv) =[I- X (0) (q, iv) X (0) (q, iv) , 

x (c) (q, iv) =\i + x (0) (q, iv) u {c) 1 x (0) (q, iv), 


(B5) 


where x^’^ (q, iv n ), X*' 0 '* (q* iv n ) and U^ s,c i are operated as 9 x 9 matrices (the upper or lower two indices are viewed as one 
number) with elements of the matrix U^ s,c i to be 


U, h — I2 — h — h 
Jh ,1 1 = ^2 / ^3 = 4 
Jh — (3 7 ^ ^2 — (4 

V, ii = (4 / I 3 = h 


(B 6 ) 


f U, h — I2 — h — h 

_ J “ 1 Jh, h = I2 ^ (B 7 ) 

lsli I Jh, h = h 7 ^ h = h 

l 2 Jh ~ V, h = h 7^ (3 = h 

The Feyman’s diagram of RPA is shown in Fig. 15 For repulsive Hubbard-interactions, the spin susceptibility is enhanced and the 
charge susceptibility is suppressed. Note that there is a critical interaction strength U c which depends on the ratio Jh/U. When 
U > U c , the denominator matrix I — x' 0> (q, iv) U i:S> in Eq.(B5) will have zero eigenvalues for some q and the renormalized 
spin susceptibility diverges there, which invalidates the RPA treatment. This divergence of spin susceptibility for U > U c implies 
magnetic order. When U < U c , the short-ranged spin or charge fluctuations would mediate Cooper pairing in the system. 

Let’s consider a Cooper pair with momentum/orbital (k'f, —k's), which could be scattered to (kp, — k</) by exchanging charge 
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or spin fluctuations. In the RPA level. The effective interaction induced by this process is as follow, 

v ?f P f A = ^ E r^(fc,fc04(k)4(-k)c s (-k0c t (k0, 

pqst,kk' 


(B8) 


We consider the three processes in Fig. 16 which contribute to the effective vertex where (a) represents the bare 

interaction vertex and (b) (c) represent the two second order perturbation processes during which spin or charge fluctuations are 
exchanged between a cooper pair. In the singlet channel, the effective vertex rfj (k, k') is given as follow. 


n / y(c) _i_ 3{/( s ) \ pt i 

r^(fc,fc') = ( 4 j +J L 


3 U {a) x {s) (Jfc - k') U W - U (c) x {c) (k - k') U {c 

1 r 


pt 


4 L 


3(7 (s) x (s) (k + k') C/ (s) - (7 (c) x (c) (k + k') U {c) 


+ 

ps 

1 

qt 


while in the triplet channel, it is 


r = 


[/(c) _ uM \ pt 1 


4 L 


U (s) x (s) (,k - k') U {s) + U {c) x (c) (k - k') U (c 


pt 


qs 


4 L 


(7 (S) X (S) (k + k') U (s) + U {c) x (c) (,k + k') t/ (c 


ps 
J qt 


(B9) 


(BIO) 


Notice that the vertex T J ^(k, k') has been symmetrized for the singlet case and anti-symmetrized for the triplet case. Generally 
we neglect the frequency-dependence of F and replace it by V x "}(k, k') « V 1 "} (k. k', 0). Usually, only the real part of F is 
adopted 25,26 

Considering only intra-band pairings, we obtain the following effective pairing interaction on the FSs, 


Veff = jj E ^(k,k04(k) C t(-k)c,(-k0c,(k'), (Bll) 

o:/3, kk' 

where a/p = 1, • • • , 3 are band indices and V“^(k, k') is 

U a/3 (k, k') = ^ r^(k,k',0)^'*(k)^>*(-k)ef(-k')ef(k'). (B12) 

pqst, kk' 


From the effective pairing interaction (Bll), one can obtain the following linearized gap equation 25,26 to determine the T c and 
the leading pairing symmetry of the system, 


1 

(2tt) 3 



V a0 {k,k') 

4 ( k ') 


A / g(k / ) = AA Q (k). 


(B13) 


This equation can be looked upon as an eigenvalue problem, where the normalized eigenvector A a (k) represents the relative 
gap function on the a —th FS patches near T c , and the eigenvalue A determines T c via T c =cut off energy e -1 ' A , where the cut 
off energy is of order band-width. The leading pairing symmetry is determined by the largest eigenvalue A of Eq.(B13). 


2. With on-site SOC 


When the on-site SOC described by Eq.(18) turns on, the free electron part of our model becomes, 


-ffband — #TB + HsOC- 


(B14) 


Let’s rewrite the interaction part in the following formalism. 


Hint 


E 




ipisO(;<TiCr2 


(B15) 
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Here Ug^ can be defined to satisfy the following symmetry relation. 


t _ rrP^f — 

u etf ~ u ea — u ez 

TT^t _ TTH v i — T/t 1 " 

U 9£l ~ U 6 £t — V 9£ > 


where the U, V tensors are obtained as 


uj 1 ! 2 = 

L 3 I 4 


v , h } 2 = 

I3I4 


lh ^ bL M=h^h = h 

h = h = h = h 
f,h=l 2 ^h = h 
~§~,h = ^3 7^ h = h 
~§~,h = h ^ I3 = h 


(B16) 


(B17) 


Due to the S z -conservation of the on-site SOC, let’s define the following six susceptibilities, 

Xh,i 2 ,i 3 ,u (q,-r) = J2 ( T - r 4 1 t( kl ’ r ) c '2t( k i + q I ' r ) c ;tt( k2 + q,°)c 4 t( k 2 ,o)) 

ki,k 2 

Xh,i 2 ,h,U (B’ 1 ") = ^ ( Tr 4 1 t( kl ’ r ) c '2t( kl +q^ T ) c ;t4.( k 2 + q,°)c 44 .( k 2 , 0 )) 

ki,k 2 

x?i,i 2 ,i 3 ,i4 (q> r ) = ^ ( T T c L( k i» T ) c <24.( k i + q> r ) c ttt( k2 + q> 0 ) c ut( k 2 , 0 )) 

ki,k 2 

h,h,h,h (q. r ) = ^ H ( r r4 l4 (ki,T)c; 2 4.(k 1 + q,r)c^(k 2 + q, 0)c ul (k 2 , 

ki,k 2 

h,h,h,h (H’ r ) = Jf Y1 (r T c| lt ( k i,r)c ;24 .( k i +q,r)c ^( k 2 + q,°)ci 4t ( k 2 , 

ki,k 2 

h,i 2 ,h,h (H’ r ) - ^ 5Z ( T r4 i; ( k i, r)c i2t ( k i + q, r)c + T ( k 2 + q, 0)c w (k 2 , 


xi 


Xh 


ki,k 2 


For (U, V, Jh) = (0, 0,0), we find the following formula of the bare susceptibilities. 


Xi 


1(0) , . , 1 


° 2 ,* 3 ,t 4 (q» iuJ n) = G £ii ( k ’ iUJ rn) G® (k + q, iu] m + «*>„) 


kjicDr, 


= Ct( k )Ct*( k )4( k +<t*( k +B)- 


N 

k,o:,/3 

2 (°) („ „v, \ - 3 ( 0 ) 


M £ k+q) - M £ k) 


ZW„ + £ k ? k +q 


Xl 1} l 2 ,l 3 ,l 4 “ Xl lt l 2 j 3t i 4 ~ 0 

XfeJLs.ufa’**'") = Xj 1 1 ( 2;l3,j 4 (q» -*Wn) 


Xi 


6(0) , . , 1 


? 2 ,; 3 ,i4 (q» «*>*») = ^ H G Ui ( k ’ G tS ( k + iu} m + iw n ) 


kjicDr, 


1 ra.*ii ua n \t-P n i \t/3,*n , -, nF ( e k+q) n -F’( £ k) 
= Tr z, 64 ( k )^if( k )Ct( k+ q)Ct( k +q)———— j — 




N 

k,a:,/3 

_ 6 ( 0 ),* 


ZW„ + e k ^k+q 


(B18) 


(B19) 


When interactions turn on, we obtain the following Dyson equations for these susceptibilities in the RPA level, represented 
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by the Feyman’s diagram shown in Fig. 15 

X 1 = X 1(0) - X 1(0) (4f7x 1 + 2V X 3 ) 

X 2 = -X 1(0) (4 Ux 2 + 2V X 4 ) 
x 3 = -X 4(0) (4 Ux 3 + 2VX 1 ) 

X 4 = x m - x m (4 Ux 4 + 2V X 2 ) 

X 5 = X 5(0) - X 5(0) (4 U - 2V) x 5 

X 6 = X 6(0) - X 6(0) (4 u - 2V) X 6 . (B20) 

Here U and V are operated as matrices with the upper or lower two tensor indices taken as one number. These equations are 
easily solved. 

The effective interaction vertex 1 ' p s '}^' 2 (k, k') is obtained from the summation of the three processes listed in Fig. 16, 

r££ ( k ’ k ') = 71 (k, k') + r^JT 1 (k, k') + r^T 1 (k, k'). (B2i) 

Firstly, the contribution from (a) is, 

r^ ?<71 (k,k ’) = {uy q %\. (B 22 ) 

Then the contribution from (b) and fc) for ft and is, 

rg Mt (k, k') = - [ 16 Ux 2 (k - k ')U + 4Vx 3 (k -k')V + SUx 1 (k - k') V + 8V X 4 (k - k') U]* 

r® P9t (k,k') = - [lGUx 1 (k - k') U + 8 Ux 2 (k - k') V + 8V X 3 (k - k ')U + 4V X 4 (k - k') V]* 

rg P9t (k, k') = 4 [(V - 2u) X 6 (k + k') (V - 2u )] p q s t . 

rg. P9t (k, k') = [16 Ux 1 (k + k')U + 8 Ux 2 (k + k')V + 8H X 3 ( k + k')U + 4V X 4 (k + k') V]^ s . (B23) 

The T p s 1°^ (k, k') for other spin configurations can be obtained by the relation T^ 9 ^ 1 (k, k') = [r^ 9 ^ 1 (—k, —k')] + from time- 
reversal symmetry. 

Finally, after the effective pairing potential (k, k') is obtained from 1 (k, k') via the relation (B12), we should solve 

the linearized gap equation (B13) for each spin configuration ■ 



